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Abstract 

Bound states of two weakly interactive massive particles are stud- 
ied. It is assumed that the WIMPonium is formed due to the gravi- 
tational interaction, since the weak interaction can sometimes be re- 
pulsive. The lifetimes of the spontaneous emission of gravitational 
radiation and of the WIMPs annihilation into a pair of gravitons are 
computed, and are shown to be many orders of magnitude larger than 
the age of the universe. 



1 Introduction 



There is accumulated evidence both from astrophysics and cosmology that 
about 1/4 of the energy budget of the universe consists of so called dark 
matter, namely a component which is non-relativistic and does not feel the 
electromagnetic nor the strong interaction. For a review on dark matter see 
e.g. [1]. Although the list of possible dark matter candidates is long (for 
a nice list see e.g. [2]), it is fair to say that the most popular dark matter 
particle is the LSP in supersymmetric models with R-parity conservation [3]. 
The superpartners that have the right properties for playing the role of cold 
dark matter in the universe are the axino, the gravitino and the lightest neu- 
tralino. By far the most discussed case in the literature is the case of the 
neutralino (see the classical review |1]), probably because of the prospects 
of possible detection. Up to now the nature of the cold dark matter par- 
ticle is unknown. Several indirect detection experiments (see e.g. [5]) are 
looking for annihilation products of dark matter in astrophysical contexts. 
Recently there has been some interest in annihilation of WIMP dark mat- 
ter via intermediate long-lived WIMPonium bound states, which can lead to 
potentially very substantial enhancements in the dark matter annihilation 
rate [6j. In these works the authors have considered the case in which the 
WIMPonium bound states are formed due to the weak interaction. This may 
sound as the only possibility, since in particle physics gravity is usually ig- 
nored, and WIMPs do not interact neither strongly nor electromagnetically. 
However, it is known from field theory |7] that if an interaction is mediated 
by a spin one gauge boson it can be either attractive or repulsive, and more 
presicely the force between two identical particles (e.g. electron-electron or 
positron-positron) is repulsive, while the force between non-identical parti- 
cles (e.g. electron-positron) is attractive. In the present work we wish to 
study the case in which WIMPonium is formed because of the gravitational 
interaction. It is interesting in this case that everything can be computed 
explicitly, and there are certain expressions for probability rates, cross sec- 
tions, lifetimes etc. As a limited case of applicability, one could have in mind 
the neutralino case (majorana fermions). 

A remark is in order here. We shall be considering a framework in which 
the weak interaction cannot bind WIMPs together to form WIMPonium. 
However, the weak interaction can still mediate transitions and annihilations 
that most likely dominate over the gravitational ones computed in the next 
sections of this article. Therefore, in reality we shall be studying of 
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academic interest in which the gravitational interactions for the aforemen- 
tioned processes dominate for a WIMP with no weak interactions at all, and 
not a generic WIMP or neutralino. 

Our work is organized as follows. The present article consists of three 
sections. After this introductory section, we present the analysis and our 
results in section 2, and finally we conclude in the last section. 



2 Analysis 

This is the main section of the present article. In the first subsection we 
recall the quantum mechanical problem for the 1/r potential case, then we 
compute the transition rate between two arbitrary bound states A and B, 
and finally in the third subsection we compute the WIMPonium lifetime due 
to WIMPs annihilation into a pair of gravitons. 



2.1 Energy levels and wavefunctions 

In the non-relativistic limit the two WIMPs (particle-antiparticle) system is 
described by the Hamiltonian 

H^^ + ^ + V(r) (1) 

where M is the WIMP mass, pi, 1 = 1,2 are the momentum operators, and 
V{r) is the gravitational potential given by 

V{r) = -G^ = ^ (2) 

with r being the distance between the particle and its antiparticle. Ignoring 
the center-of-mass motion, which is that of a free particle with mass 2M, the 
dynamics of the two-body system is described by the Hamiltonian 

h = f^ + Vir) (3) 

where — M/2 is the reduced mass of the system, and p is the momentum 
operator conjugate to the position operator r. Because of the spherical sym- 
metry it is convenient to use the spherical coordinate system, in which case 
one has to solve a one-dimensional problem with respect to r 

y"{r) + 2i,{E-V,ff{r)) = (4) 
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where the effective potential is given by 

Ven - y{r) + (5) 

with / being the angular momentum quantum number. Finally, the wave- 
function is given by 

^nim(r,e,<i>)^Yr(e,<i>f-^ (6) 

with Yj^ the spherical harmonics, and with the normalization conditions 

1 = / dn\Yi"'\'' (7) 
1 



POO 

/ dr\y\' (8) 

^0 



while the allowed energy levels are given by 

En^^ (9) 

where n is the principal quantum number, and Ei is the energy corresponding 
to the ground state. Later on we shall make use of Ei,ao, and therefore we 
report here their values 

ao = ^ (10) 
nag 

El = (11) 



or equivalently 



ao = ^{^yiri (12) 



El = -T ] M (13) 

where rupi ~ 10^^ GeV is the Planck mass, and Ipi ~ 10~^^ cm is the Planck 
length. 
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2.2 Spontaneous emission of gravitational radiation in 
the quadrupole approximation 

In the case of spontaneous emission, a WIMPonium state A makes a tran- 
sition to a state B in the absence of any incident gravitational wave. In 
this subsection we shall derive the lifetime of the state A by exploiting the 
similarities with the electromagnetic theory. Let us first recall here that in 
atomic transitions the dominant contribution (whenever is allowed by the 
selection rules) in the spontaneous emission of electromagnetic radiation is 
the transition in the electric dipole approximation. In the quantum theory 
of radiation one obtains the following expression for the integrated transition 
probability between two atomic states A and B [8] 



A->B 



or 

a; , _ 12 



xbaI' (14) 



i- A^>B = \Pba\ (15) 

OTT 

where e is the electron charge, u is the energy of the photon emitted, xba = 
{B\x\A) is the matrix element of the position operator between the states 
A and B, and we have defined p = ex. On the other hand, in the classical 
electromagnetic theory the average power for the dipole radiation is given 
by|g 

^ = ^7 

i/vr 

where Pq is the amplitude of the oscillatory electric dipole moment of the 
radiating source, p{t) = pocos{ut). Remarkably enough, one can use the 
classical result to obtain the formula for the transition rate rA->B given 
above, avoiding the full quantum mechanical computation. Within the time 
interval r = T~^, with r being the lifetime of the state A, one photon with 
energy u is emitted, and therefore we can write 

P = - = uT (17) 

T 

Therefore, the transition rate F is determined and given by 

r = ^ = ^|^r (18) 
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which means that from the classical result, the quantum mechanical result 
is obtained provided that we replace po/2 by the matrix element Pba- This 
is not strange at all, if we recall that in the classical treatment the radiat- 
ing source is characterized by an oscillatory dipole moment pit) = pocos{(x)t), 
while in the quantum theory of radiation the interaction Hamiltonian is writ- 
ten in the form Hi{t) = Hjexp{±iut) . Therefore, if we write the cos{ujt) as a 
sum of two exponentials, cos{ut) = (l/2)(exp(ici;t) + exp(— zwt)), we recover 
the substitutional rule po/2 —> pba- 

Now let us move to the gravitational radiation case. It is known [10] that 
in the classical theory of general relativity the emission of radiation due to 
electric or magnetic dipole vanishes. The first non-vanishing contribution 
comes from the quadrupole emission mechanism. So, we need the classical 
result for the power for gravitational radiation in the quadrupole approxima- 
tion, and then we can easily obtain the corresponding quantum mechanical 
result. In the classical electromagnetic theory, the power for radiation in the 
quadrupole approximation is given by |TT] 



p-m='-i^j:\q.,\' (19) 

where kc = l/(47reo) is the constant entering into Coulomb's law for the 
electric force between two charges, and the quadrupole tensor Qa/3 is defined 
as 

QaP = j d^Xp{x){'iXaXp -Ixl'^ 5ap) (20) 

with p{x) being the electric charge density with the property 

d^xp{x) = Qtot (21) 

where Qtot is the total electric charge in space. The classical result for 
gravitational radiation in the quadrupole radiation looks like the one for 
electromagnetic radiation 



, \^al3\ 

a,/3 



^ 360 ^ 



/a/^r (22) 



where now k^. is replaced by the Newton's constant G, the extra numeri- 
cal factor of 4 is due to tensor calculus pi]], and /q,^ is the gravitational 
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counterpart of Qa/3- Now we have the mass density with the property 

J (fxpuix) ^ Mtot (23) 
with Mtot being the total mass in space, and the mass quadrupole tensor 

lafi = J d^XpM{x){?>XaXp - 5ap) (24) 

and for the WIMPonium an order-of-magnitude estimation is 

/a/3 - Mai (25) 

Finally, by applying the rule /q^/2 — > we obtain the transition rate 
between two WIMPonium states A and B in the quadrupole approximation 

r^:,, = ^ (26) 



45 

a,/3 



where the matrix element reads 



I^A = (27) 

Taking into account that G — l/m^i, ui ~ and {B\Iap\A) ~ Moq, we 
finally obtain the formula 



rA->B-(^) t,, (28) 

where tpi ~ 10~^^ sec is the Planck time. Assuming a WIMP mass M — 
10 TeV we find for the lifetime 

TA->B ^ (29) 

with tQ being the age of the universe, and we have used the numerical values 

to ~ 10^°7/r (30) 
yr ~ 3 X lO^sec (31) 



15 
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2.3 WIMPonium ground state annihilation into a pair 
of gravitons 

We assume that the WIMP particles responsible for the cold dark matter in 
the universe are scalar fields. The available data do not require that WIMPs 
must be fermions, and as a matter of fact there are already in the literature 
scalar field candidates for playing the role of dark matter. We can mention 
at least the case of the axion [12j , which has been introduced in order to solve 
the strong CP problem via the Peccei-Quinn mechanism [13], and the case 
of branons that are related to the fiuctuations of a D-brane [H] on which we 
are confined according to the brane- world idea [.15J. We choose to work with 
scalar fields to avoind extra complications from spinorial fields, but our final 
results should also hold for fermionic WIMPs, since we eventually take the 
non-relativistic limit. 

First let us recall how one computes the annihilation rate of the positro- 
nium ground state into a pair of photons. Ignoring for one moment that we 
have a bound state, one needs to compute the total unpolarized cross section 
of the process, e~e"'" — >■ 77, using standard field-theoretic techniques, then 
take the non-relativistic limit, and finally take into account that what we 
have is really a bound state. For that we also need the ground state wave- 
function, which has already been determined by solving the Schoedinger's 
equation. All in all, the inverse of the positronium lifetime (ground state) is 
given by |8j 

Tin = 1,'S^ 27) = AaZr'v\^is{x = 0)\' (32) 

where the factor of 4 is due to the spin, namely we have taken into account 
that in the positronium ground state the single state does the whole job. 

Now let us move to the graviton- WIMPs system. According to general 
relativity, if we have a canonical scalar field coupled to gravity, then the 
total action describing the model contains the Einstein-Hilbert term for the 
gravity part plus the appropriate term for the scalar field 

S = j d'x^g A + 1(90)2 - V{<P)^ (33) 

where R is the Ricci scalar, g is the determinant of the metric, k = VSirG is 
the strength of the gravitational interaction, and V{(j)) is the self-interacting 
potential for the scalar field. For the discussion to follow we shall consider 
the case of a massive scalar field, V{(j)) = {l/2)M'^(f)'^ . If we write the metric 
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as 9fj.u = Vi^u + t^h^i,, it possible to obtain a usual field theory in Minskowski 
spacetime, the "scalar gravito dynamics" , which describes a massless spin two 
particle (graviton), a massive spin zero particle (WIMP) interacting with each 
other. The total Lagrangian takes the form ^16j 



S = j (fxC (34) 
C = C\h^,]+C\(P\-'^W''T^, + C[hhh] + ... (35) 

where the first term is for free gravitons, the second term is for free WIMPs, 
the third term is responsible for the interaction between the gravitons and 
WIMPs, and the rest of the Lagrangian describes the interaction of gravitons 
among themselves. Here T^,y is the energy-momentum tensor for a scalar field 

T,, = d.cpd^ - Q(90)2 - V{<P)^ g,, (36) 

Given this interaction term one can derive the Feynman rules for the in- 
teraction vertices of this model [11]. Regarding the free part of the total 
Lagrangian, one needs to quantize the free scalar field, as well as the free 
graviton. The scalar field quantization is well established in many standard 
books on field theory. As for the graviton case, the treatment should be 
similar to quantum electrodynamics. There are a few differences due to the 
fact that the photon field carries a single Lorentz index (A^, spin one boson, 
polarization vector e^^^ 2(^))' while the graviton carries two Lorentz indices 
{h^^ , spin two boson, polarization vector e^=i2(^))- Finally, the Feynman 
propagator in momentum space for a free graviton reads (in the commonly 
used De Bonder gauge) [T7] 

Dr'ik) = ^—-^ (37) 

In a similar way to the electromagnetic case and the positronium, we are now 
interested in the process, 00 — )■ GG, with G being the graviton, and we first 
need to obtain the total unpolarized cross section at tree level and in the 
non-relativistic limit. Then we can obtain the lifetime of the WIMPonium 
ground state taking into account the corresponding wavefunction. For the 
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computation of the total cross section for the process 1 + 2— i'3 + 4we could 
start from the standard expression for the differential cross section in the 
center-of-mass system given by [8J 

where s is the Mandelstam's variable, Pi and pf are the initial and final mo- 
menta respectively, and M. is the Feynman amplitude computed by applying 
the Feynman rules for a given model. Then we could apply Feynman's rules 
to compute the invariant amplitude Ai, and finally we could integrate over 
the solid angle to obtain the total cross section. However, we do not really 
have to go into this, since we are only interested in the non-relativistic limit, 
M ^ K, with K being the kinetic energy, in which the energy of the WIMP 
is ~ M, and the WIMP momentum is p ~ Mv. In this limit we can see 
from eq. ( l39l) that the cross section goes like a ~ 1/v, and from second order 
perturbation theory and dimensional analysis the cross section should be 

a ~ (40) 

V 

where the precise numerical prefactor turns out to be 27r [18]. Therefore, 
in total the WIMPonium lifetime is given by (since now there is no extra 
complication related to the WIMP spin) 

''"wiMPonium = rwiMPonium = (T N Rv]"^ {0)\'^ = 2nG'^ M'^ 3 = - ( ) M 



Trag 4 \mpij 



(41) 



or 



TWIMPonium — 4 i^'J^ j tpl (42) 

If the WIMP mass is M = 10 TeV, then the last formula for the WIMPonium 
lifetime gives 

TwiMPonium — 10^°^ to (43) 

Therefore, according to our numerical results we conclude that the lifetimes 
(both for WIMP annihilation into a pair of gravitons and for transitions 
between bound states in the quadrupole approximation) are many orders of 
magnitude (in fact astronomically large!) larger than the age of the universe. 
Since we have obtained these results in the non-relativistic limit, we expect 
the same result to be valid for fermionic dark matter particles, although our 
discussion has been focused on scalar dark matter particles. 
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3 Conlusions 



In the present work we have considered bound states formed due to the 
gravitational interaction between two identical weakly interactive massive 
particles. These particles are the ones responsible for the cold dark matter 
in the universe, and therefore they no not have neither strong neither elec- 
tromagnetic interactions. Furthermore, if the particle is identical to its own 
antiparticle we expect from general field theoretic arguments that the weak 
interaction between the particle and its antiparticle should be repulsive. In 
this bound state can be formed due to the gravitational interaction 

only, which is a universally attractive interaction. We have assumed that 
WIMPs are scalar particles for simplicity, in which case we can easily write 
down the Lagrangian for the WIMP-graviton system, and determine the 
Feynman rules for computing certain processes in quantum gravity. We do 
not worry about the non-renormalizability of quantum gravity, since in this 
work we are only interested in processes at tree level. In the non-relativistic 
limit the allowed energy levels and the corresponding wavefunctions of WIM- 
Poniums can be computed exactly, and are known from the familiar Hydrogen 
atom case. We have computed the lifetime of a bound state due to sponta- 
neous emission of gravitational radiation in the quadrupole approximation, 
as well as the lifetime of the ground state of WIMPonium due to the WIMPs 
annihilation to a pair of gravitons. We have determined numerical values for 
the lifetimes, which are shown to be many orders of magnitude larger than 
the age of the universe. 
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